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SYNOPSIS 

Aerial photographic methods are now widely used in topographic surveying, 
and in recent years efforts have been made to place them under mathematical 
control by analytical procedures of computation. In this field Earl Church,? 
Assoc. M. ASCE, and P. H. Underwood,? M. ASCE, are the leaders. The 
purpose of this paper is to present an analytical solution based on the direction- 
cosine principles of space geometry. The underlying mathematical principles 
are first developed and then applied to aerial survey problems. The methods 
of analysis and adjustment and the determination of the effect of the atmos- 
pheric refraction are new. 


Data UsED 


The data immediately available in an aerial survey are the projections of 
the ground points on the negative plate, as transmitted through the objective 


lens; but the angles at the front nodal point, subtended by the separations of 


the objects at the instant of exposure, will be used in the following analysis. 
Furthermore, the positions of three or more ground points are also required, 
and these will be provided by the ground survey. The space angles at the 


front nodal point are assumed equal to their corresponding angles at the rear 


nodal point. The latter may be determined if the perpendicular distance from 
the rear nodal point to the plane of the negative plate, and the coordinates of 


_ the base of this perpendicular with respect to the fiducial axes on the negative 


plate, can be determined. An analysis of the geometrical relation between the 
rear nodal point and the fiducial axes is a part of this paper. The effect of 
atmospheric refraction is also analyzed; but the measurement of the coordinates 


of points in the negative by a comparator or other special instrument, and the 
ETS a et aes NE I LEE EE ares 


Nore.—Written comments are invited for publication; the last discussion should be submitted by 
September 1, 1950. 

1 Prof, of Surveying and Geodesy, Tangshan Eng. College, Tangshan, North China. 

2 ‘Manual of Photogrammetry,” Pitman Pub. Corp., New York, N. Y., and Chicago, IIl., 1944, p. 536. 

3‘Space Resection Problems in Photogrammetry,’ by P. H. Underwood, Transactions, ASCE, Vol. 


112, 1947, p. 921. 
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calibration of the camera to determine values of distortion, are not considered 
in the present paper. 


MATHEMATICAL PRINCIPLES 


For convenience of reference, a summary of the mathematical principles 
directly involved in the following analysis will be presented first. Some well- 
known formulas are given without proof; but, for those not commonly known, 
the method of derivation will be outlined. . 

The method of adjustment described consists largely of the determination 
of corrections to be applied to computed or observed values. When the true 
values are connected by one or several conditional equations, the approximate 
relations between the corrections may be obtained by differentiation of such 
conditional equations, since the corrections to be applied are quite small in 
comparison with the values to be corrected. The resulting formulas will all 
be linear equations containing the corrections as unknowns and the differentials 
as exponents whose values may be found from the known values of the quan- 
_ tities. This general method will be used to derive many of the required 

equations. 

Reference Axes.—The normal line through point O on the sphere OSWNHE, 
Fig. 1, is adopted as the Z-axis; the north-south line through O, in the tan- 
gential plane of the sphere at point O, is the Y-axis; and the east-west line is 
the X-axis. The positive directions are all indicated by arrowheads as shown. 
The center of the sphere is designated point C, and its radius by R. 


Normal Distance of a Point to the Sphere—Let the coordinates of point P 
be X,, Y,, and Z,, and let a normal line from point P on the sphere be drawn 
to point P; as shown in Fig. 1. Then ‘ 


PP =(X%+ ¥%,4 (2+ RY} — BR 


X*, + VY, _ (X*» + Y?,)? 


Jet 3G, +B) 1 ae PRT ee a 


which may be used to determine the elevation of a point defined by coordinate — 


| distances. 
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Length of the Arc OP;.—The curved distance, OP,, on the sphere is given by 


DPD. ~1 (Xx?) aie ¥35)3 Saks 


P,; = R sin ais aunel PRLS AON gat HIRE CG NS 2 
f : REDD (2) 
Eq. 2 may be used to determine the length of line OP at sea level. 
Equation of a Normal Line——The equation of the radial line PC is 
IOS Caine Lie a? (3 
Yo epee a ee ) 


which may be used to determine the equation of the plumb line from a point 
having known coordinates. 

Coordinates of a Point from Its Elevation and Distance-—Let the azimuth of 
line OP, measured from the north point in ¢lockwise direction be a; its curved 
distance OP; be s; and its elevation P;P beh. Then 


X, = (R+h) sin = sin a Fut, ab ies As eae (4a) 

Vy = (B+ A) sin 5 cosa... iy. sees e vs (4B) 
and 

Yume KR A) cos BR a tee ok ee ee (4c) 


Eqs. 4 may be used to determine the coordinates of points located by a ground 
survey. 
Azimuth Angle of a Line-—Let N’ be a point on the Y-axis and let the ver- 
tical angle of the line CN’ be ¢ as shown in Fig. 1. Then the coordinates of N’ 
are 0, R cot ¢, and 0; and the equation of the plane containing points N’, C, 
and P is 


aX Y. Z 1 

0 Recot¢ 0 te es 

; 6 _R 1 mat Qe seen eee ee (5) 
Xo» ie Zp 1 


Let P, be another point with coordinates X1, Yi, and Z;. The equation of 
the plane containing points C, P, and P, is found to be 


xX Y x 1 
0 0 —R La 

os ()) eh oe ae 6 
Sete ra eae (6) 
Xy Yi Zi 1 


_ Eq. 5 can be used to determine the meridian plane through a given point; 
Eq. 6 gives the vertical plane through two given points; and the angle between 
the two planes is the azimuth of line PP, which may be determined from the 
equations of the two planes by the usual method. 
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Space Angles at a Point—The coordinates OS., OS,, and §.58 of point 8, 
Fig. 2, with respect to the rectangular axes are X;, Y,, and Z,, respectively. 


Similarly the coordinates of points P, and P2 in the (X Y)-plane are denoted by 
Xi, Y1, X2, and Y2, respectively. Angle P:SP2 may be determined from these 
coordinates thus: 


S.Pi = V (Xa — Xe)? + (V1 — Ya)? ose ese ees (7a) 
SiP2 = V(X2 — Xe)? + (V2 — Y,)?. «0 eee eee (7b) 
PyPs = -V (Xa — Aah (We = Yate lass eee (7c) 
BPi NIG ASLP ee ee (7d) 
SPs V2, (Sib) (7e) 
and ase, ea eee 
cos Z P,SP2 = (BEDE CE ARES)! ARIS) eee (7f) 


2 SP, SP. 


Eqs. 7 can be used to determine the space angle at an exposure station. Lines 
SS. and SS, are perpendicular to axes OX and OY, respectively. 

Coordinates in Terms of Axial Distances——Let points I,, Iz, and Is, Fig. 2, 
be points on the X-axis forming equal space angles I,SI2: and I.SI; to be de- 
noted by 6; and let the angle 8,SI; be denoted by Z,. Then, from the relation, 


SI, _ [ls _ cos(@—#.) . 
SI, ~ TL ee LEB.) it may be found that 


tan E, = dates GOEL ie eee Yi sae ee (8) 


II; 
Since S,I, = SS, tan F, and Sals + II; = SS. tan (6 + E,), it may be found 
that , 
Son ih 
Sees fan OPE.) otan Be (9) 
and 
nN II; tan E, 
8,1. = tan (6 + E,) — tan A, © 0 0 6s 6 06 00 be ee (10) 
From Fig. 2, it is evident that 
X,.= Ol: — 8203.02. it ee (11) 
and 
Sih era mn Te 
Z, = V(SS;)? — ¥3, = V(§8,)? — A he ie ee ee (12) 


to the fiducial axes. 


= 
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Equation of Straight Line—The coordinates of any point P, will be denoted 

by Xn, Yn, and Z, and their small corrections by 2m, yn, and Z,. Direction 

cosines of a line in the directions of X, Y, and Z will be denoted by L, M, and 

N and their small corrections by J, m, and n, respectively. The equation of 
line PP: is ' 

X—-Xi Y-¥Y, 24-42, 


[oo apres Soles pee crag ee (13) 
in which 
Pe 2 = Sv eine Me ee (14a) 
Wa? = A Eee Teens MOLT re (140) 
and 
= Neate adi An es (14e) 
In Eqs. 14, 
PW (kg ack (Yar Val (Ze = Zi een: (15) 


If small corrections x1, y1, and 2; are added to the coordinates of point Pi, 
it can be found that 


b=2[-(-DatlLMntlNal Oe en (16a) 
m=2(LMx-(1-M)y+MNaj eae ee (168) 
“and 
=*(LNa+MNy— (1-4) a] Free ieee bt (16c) 
_ Eq. 16c¢ can be written 
Lli+Mm 
Za ss shan eh vee (17) 


Distance from a Point to a Line—The perpendicular distance D, from any 
point P, to the line defined by Eq. 13 is given by 


D2, = (Xi a Xe)? + (V1 ip RS 
PAs. = Zp) ah (Xe = Xs) 
+ M(¥, — Yn) + N (2: — Z,)}-. . 18) 


Angle Between Two Intersecting Lines.—Let the 
direction cosines of two lines SP: and SP2, Fig. 3, be 
denoted by Li, Mi, Ni, L2, M2, and Nz, respectively. 
‘The angle Aiz included by the two lines (P:S and 
SP2, Fig. 3) is given by 

cos Ave = 1, L22+ Mi M2.+NiN2 Sihsttel Shenae leva eveus (19a) 


The correction ai (in radians) of angle Ai2 due to small corrections of the 
direction cosines is given by 


G12 = — cosec Ayz (Lol; + Lile + Mzmi + Mi mz + Nam + Ni me). . (190) 


a 
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. d 
If the corrections of the direction cosines are due to the corrections 2, y, 


and z of the coordinates of point S alone, it may be found that 


pe eee Ai {C(ri Li + re Le) — (71 Lo + 72 Li) cos Aiz] x 


T1172 


ae [(r1 My, + r2 M2) a (ry M.+ To M) cos Ain] y 
=F [(ri Ni + r2 Ne) at (r1 No + r2 Ni) cos A12] 2}. . (20) 


in which r; and rz are the original lengths of the lines SP; and SP2, respectively. 
Distance Between Two Lines——The perpendicular distance between two 


‘ : lines. eee eek ae es Fp aaee.< ee le 
nonintersecting lines, AED I pee Ta an tee 
wit ee is given by ° 

X,— X2 Y,—-—Y2 Zi — Ze 
Ty M, Ni 
Tz M2 N2 


it STM, Ny = ata NG! PAL Na La Ni + Ma Le 


The condition of intersection of the two lines is controlled by 


Xi — X2 Y,-—Yf2 Zi — Ze 
Ty mM, Mi 
Lp M, N2 


Small corrections 21, 22, yi, Y2, 21, 22, li, le, m1, M2, 21, and nz to be added to 
the coordinates and the direction cosines of two nearly intersecting lines to 
make them intersect are connected by 


X,— Xe Y,- Yo Z1— 22 X, — Xo Y,;- Yo 21 — 22 
Ty M, Mi = Ty My, ny 
Le M, No Le Me Ne 
Xi — X2 Yi-Y2 Z1—Ze %i—% Y¥i—Y2 21-22 
+ Ly ™m, Ni + ly My, N, = 0), . (22) 
Le M2 N2 ls M, N2 


Kq. 22 was obtained by neglecting small terms involving powers of the correc- 
tions, which powers are higher than the first. 

The coordinates Xm, Ym, and Zm of the middle point P,, of the line perpen- 
dicular to the foregoing two lines may be found to be 4 


Xm a 3 (X1 + X2) 


Xi 53 XG Y, => Yat Z1- Zs 
2L; Le IiM:+ L2M, Ii No+ L2N, 
M,N2—M.N, I2N,—1L,N2 Ii M2 — Le My 


200 - MP + aM — Lh N2)? + (Li M2 — Le M,)?) Gay 


} 
= | 
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Ym = 4(¥1+ Y¥2) 
X, — Xo Yi-Y2 Z1—Z2 


LI, Mz + Le M, 2MiM, | MiN2+M.N, 
M,N2— M.M, L2N, — 1, N2 Li M2 — Le My, 


7 20s Na = Ma Ma (i Ni — Ly 2) + (Ly Ma — Le My Py 2) 


and 
Zm = % (Zi + Z2) 
X;— Xo Yi-— Y/Y. Zi — Ze 
[Li Ne+ L2Ni M,N.+ M2N, 2NivNz2 | 
M,N2— M.N, L2N, — 11 Ne Ii M2 — Le My, 


5 2((Mi N2 — Mz Ni)? + (L2 Ni — Li No)? + (Ly M2 — La My)? . (23c) 


which may be used to determine the most probable point of intersection of two 
observed lines that should meet at a point if there were no errors of observation. 

Three Concurrent Lines.—Let the space angles between the three concurrent 
lines SPi, SP2, and SPs, as shown in Fig. 3, be denoted by Ais, Ais, and Ags, 
respectively; and let H, denote the dihedral angle between the two planes 
SPiP2 and SP:P;. By the law of cosines in spherical trigonometry, 


cos A23 — cos Ai2 cos A413 


neo fi, = sin Ais sin Ais 
If H; remains unchanged, but Ai. and A1i3 are corrected by ai2 = — ptan Aj2 
and a1; = — p tan Aj;3, respectively, the’ correction to be applied to Ag; is 
given by 
cos Ai3-~—s coS Aye 
Qos = p cosec Azg (2.co8 Ass aioe A aA seit Wie seat eO) 


When the direction cosines Z;, My, Ni, L2, Me, and N2 and the space angles 
Ai3 and A23 are known, the direction cosines of line SPs; may be found by 
solving the equations: 


cos Ayi3 = 2,23; + MiM3;+ Ni N3 re a ai etitl ecg ee (26a) 
cos Aes = Le L3 + MeM3+N2N3 Vai ok Soto: (26b) 

and 
LD; + M?, + N?; a1 We nie AALS Laie ver eae ve | (26c) 


From Eqs. 26a and 268, the magnitudes L; and M; may be found in terms of N3; 
and, by substitution in Eq. 26c, a quadratic is obtained which may be solved 
for the root N3. 

If points P; and P, are fixed in position, if point S is moved a little to a 
new position having (X + x), (Y + y), and (Z + 2) as its new coordinates, and 
f the space angles A13 and Az; are given small corrections a1; and a3, respec- 
tively, the corrections li, mi, 71, 12, m2, and m2 may be found in terms of 2, y, 
and z by Eas. 16 and 1, ms, and n; may be found in terms of 2, y, and z through 
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the solution of the equations: 


Lils + Mims + Nins = — sin Ais a13 — L3l, — Ms m — Ns m..(27ay 


Lolsy + M2m3+ Neons ae sin Ao3 23 — D3 lz — M3 m2 — N3 Nz. . (276) 


and 


which are obtained by differentiating Eqs. 26. 


Four Concurrent Lines—Let the space angles between the four lines SPy, 
SP, SPs, and SP, in Fig. 3, be denoted by A1, Aj3, Au, Aos, Aa, and Asa 
Since the dihedral angles between the planes SPiP2, SP,Ps, and SPiP, can be 
expressed in terms of the space angles by Eq. 24 and the sum or difference of 
two is equal to the third, by equating the cosine of the one to the cosine of the 
sum of the other two, it is found that 


Ky23 K 
cos Az, = cos A13 cos Aig + Sean, 
4 kc08 Ass = cos:A1s.¢oa Ais) (008 Ase 00s 415-008 2 in) aan 


sin? A 12 
in which 


Kiss = (1 + 2cos Ai2 cos A13 cos A2s — cos? Aiz — cos? Ais — cos? Aas). . (29a) 
and 
Ki. = (1 + 2 cos Ay cos Au cos Au — cos Ay ee cos? Ay. — cos? A o4)?. . (296) 


The second term on the right side of Eq. 28 will have a negative sign wh 
points P; and P, are on the opposite sides of plane SP:P2. The small correc- 
tions of the six space angles are connected by 


G34 = Ci2 Ai2 + C13 Giz + C14 Gig + C23 Gog + Cos Org... 20s (30) 


in which the differential coefficients are 


Ky (cos A1z — cos A13 cos A23) 
Cig = 2 cot Ais'cot Ag — ' 
3 ° Ki03 sin Ax sin Aza 
a Ki23 (cos Ai — cos Au cos Ax) 
Ky sin Ai sin Asa 


__ cos Ai, cos Avs + cos Ai3 cos Ang 


ae Fy BF - - (31a), 


_ Sin Ais (cos Ais — cos Ai2 cos Ax) 


C13 ; ; 
sin? Aj. sin Ag, 


_ Kies sin A13 (cos Ais — cos Ais cos A23) 
K423 sin? Ai sin Asa 


. (316) 
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sin Ais (cos Ais = cos A12 cos A2s) 


C14 = —> 
sin? Ajosin Ags 
P 
_ Kie3 sin Ais (cos Aig — cos A12 cos Axa) (31¢) 
Ki. sin? Ai sin As at c 
oe sin A23 (cos Aes — cos Aiz cos Aj4) 
sin? Ajie sin Ag, 
_ Kia sin Ags (cos A2s — cos A12 cos A13) (314) 
K 423 sin? Ais sin Aza ie 
and 
a _ sin Ax (cos A2s — cos Ai cos Ais) 
: sin? Aj. sin Ag, 
ine Ky23 sin An (cos Au — cos Ay cos As) (31 ) 
Ble 


Ky sin? Axe sin A3z4 


If the angles Axo, A13, Ais, Aes, Aes, and A3y are determined by the coordi- 
nates of points 8, Pi, P2, Ps, and Pu, the value of Ags, obtained from the triangle 
SP3P., will be identical with that determined from Eq. 28. 


If the conditional formula (Eq. 28) is satisfied, the direction cosines Ls, Ma, 
‘and N, of line SP, are determined from the values of the direction cosines in 
the other three lines by the simultaneous linear equations: 


cos Ay = L,l4+M,M,t+Ni Ns GRE Palm comand ie Cayce (32a) 

cos Ang => Le Is + M, M, + No Ni Tesla pet es Sees eheehexs (326) 
and 

cos Asa = L3 La + Mz M4 + N3 Na dels icawier alte nusllquarehe (32c) 


Angles and Sides of a Tetrahedron.—The space angles of the tetrahedron at 
the vertex S, Fig. 4, are denoted by Ais, Ais, and A2zs; and the lengths of the 
sides, SP, SP», SP3, P,P., P;Ps, and P2Ps, are denoted by D,, Dz, Ds, Dy, D,s, 
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and Dos, respectively. The relations between these angles and sides are 


Diy D8 Dt 2 Da Ds C08 Ais e on een eee (33a) 

D?*13 = D*, + D*; —2 dD, D3 cos Ais So Sr ae erase & Doe (336) 
and ; 

D3 = D*, + D?; —2 Dz D3 cos Ao acaliele elec ig Tha he Canoe (33c) 


which are usually used to determine D;, Ds, and D; from the known values of 
Ais, Ais, Ag3, Dio, Ds, and Dy;. The small corrections di, do, and ds; of D,, 
D2, and Ds, respectively, are connected by 


— Di — D2 cos Arg 
D, cos Ais ‘Se D2 Bd 0 0 celle’ © ©. -ei'u te 0X lip) esipahie 
— Di — Ds cos Ais 9 

D; cos Ai3 — Dy ae 


de 


and 
Ds; = Dz cos Ao3 ds 


dz © Ds Cos Asses at ae He eee © <a bw ee ye OH Ouro. 


Combining Eqs. 346 and 34c, 


1, = Dia Dscos Ais. Ds — D2 cos Ars 9 
~ Dy cos Ais et D3; Ds; cos Ao aa D» wa: a Sete“) ate ails! 


Geometrical Solution of a Tetrahedron.—To determine D,, D2, and Ds; geo- 
metrically, first construct a plane triangle P,;P2P; to scale and draw circular 
ares on its three sides to contain the vertex angles Ai2, Ais, and Ags, respec- 
tively, as shown in Fig. 5. Then with an estimated value of D, as radius and 
point P, as a center draw a circular arc to intersect the circles on either side of 
point P, at points B; and C; as shown. Next, with P2B: (= P3:C:) as radius, 

draw the intersecting arcs at 11, Fig. 5, using points P2 and Ps, respectively, 
as centers. If point I, falls within the third circle as shown, the assumed 
value of D, is probably too long. 

Use a revised value of D, to repeat the same operation, biked determining 
another point Iz, and, if necessary, a third point I3. Let the smooth curve 
through points I, Iz, and Is intersect the outer arc of the third circle at point I, 
Fig. 5. The distances P2I and P;I are the required values of D2 and D; 
respectively. With line PI representing the value of Dz, it is an easy matter 
to determine D;. The foregoing method can be used to determine D,, Dz, and 
D; to three decimal places. 

Analytical Solution of a Tetrahedron—The approximate value D’, deter-— 
mined by the geometrical method is used for D, in determining the value D’2 
of D2 by Eq. 33a and the value D’; of D; by Eq. 33b. Then D’; is introduced 
as the value of D; to find the value D’’: of Dz by Eq. 33c, and finally the cor- 
rection d; is found by 


Di s=.D!, cos Ais) 2 py DV SAD Gos Ags 
Berens Aties Da d= D a Os DY coe Ang ade 


D'; — D": cos Ars | 


x T0088 Au =D 


. . (36) 
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Eq. 36, solved for d,, becomes 
foe = D”, 


ey MS PUA TEM er Reece nd dO a RY 7 7 
ah D'; — D's cos Ais D's Coy Dir cos Ao Lid D'; — D', cos Ax ( ) 
DD; cos Ais — Dre D's cos Ao mae D5 D'; cos Ai => D', 

- The values obtainable from the foregoing process of calculation are 

De Dyck duis die weet eae (38a) 
{ es ’ D'; = D’, COs Ais 
Dy, = D', + Diioea Dh, Geilo ce Sa rea (38d) 
and 
JT y/, 

jay Og ra reac Ne SE ee A (38c) 


D'; cos Ai3 — D'; 


More accurate values can be found by using the improved value of D, to 
repeat the process. 
If the coordinates Xi, Xoo “Li, Xo, Yo, Za, X3, Ys, and Zz of points Pi; P., 
and Ps, and the lengths of the sides D,, D2, and D; are known, coordinates X,. 
Y, and Z of point S are found by solving the equations: 


DP CEN 1) a (ee (ona) ee al ee (39a) 
D*, = (X — X2)? + (Y — Yo)? + (Z — Ze)?.. 22... (39d) 
and 


D?; = (X — X3)? + (Y — Ys)? + (Z — Z3)?........ . (89) 
from which it follows that 


BG — X) X +2(%1— Ys) ¥+2(Zi— 2) 2 
3 = X*%, + Y%, + 24, — X% — Y%. — 2%, — D*, + D*2. . (40a) 
and 
2(Xi — Xs) X+2(¥1 — Ys) ¥+2(21—-—2;)Z 
= X?; + Y*, + VA rae XxX? ev Y?; Tis. 2’; vies D*, + D?;. . (400) 


Thus, X and Y are found in terms of Z in the following form: 


and 


in which C1, C2, ki, and k, are numerical constants. Substituting the equivalents 
‘of X and Y into one of the original equations and solving it, two values of Z 
are obtained which correspond to two positions of point 8, one above and one 
below the base P:P2Ps. 

Refraction of Air—According to the theory of atmospheric refraction, the 
total change of direction of the light ray from its position at infinity to point 
§, in the atmosphere, as shown in Fig. 6, is given approximately by 


eat (Gin V) CAC is Sates ce ee reer (42a) 
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in which C; is the angle between the tangent line of the path of light at point 
S, and its original direction; {1 is the angle between the tangent and the plumb 
line; and g; is the refractive index of air at point S:. Similarly, at another 
point Sz on the same path of light, 


Co = (g2 — 1) tan fe..... s 5 wide haat ueat ieee (42b) 

The change of direction of the curve from points 8; to S2 is equal to 
Cis = (g2 — 1) tan £2 — (g1 — 1) tam f1...........-- (42c) 
If the path of the light ray is very flat (as is actually the case), [1 and f2 may be 


F, 


Fia. 6 Fie. 7 Fie. 8 
considered as equal, and, therefore, 


Ci. = (g2 — 91) tan 1 és ere getel ly wile be 6 et Slecdn sens eee (42d) 


In that case it may be reasonably assumed that the angle between the chord 
§.S2 and the tangent at point 8: is proportional to C12. Consequently, 


Z TSiS82 => — ptan fy 


in which p is an unknown constant whose value mainly depends on the height 
of point S; above point S:. The negative sign is used to account for the fact — 
that the zenith distance of the direction of the chord 8.8; is smaller than £1, the — 
direction of the tangent TS:. Since angle TS,S: is always less than the astro- 
nomical refraction in numerical value, the value of p in minutes of arc is less’ 
than 1, if the zenith distance does not exceed 75°. 


PROBLEMS IN AERIAL SURVEY 


Camera Constants——The relation between the rear nodal point and the 
fiducial axes on the negative plate is essentially as shown in Fig. 2. If the — 
images I, Iz, and I; of three objects which are at equal angles from each other - 
appear on both fiducial axes, the actual angles 6 at the exposure station between 
the real objects are measured with a theodolite, the lengths Ol, Ol, and OIs 
along the fiducial axes on the negative plate can be determined accurately 
(corrected for distortion if necessary), and the three coordinates of the rear 
nodal point may be calculated by Eqs. 8 through 12. On account of the errors 
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of measurement and distortion, the values determined from different sets of 
data may be slightly different. 

When the coordinates of the rear nodal point are exactly known, the space 
angle at the node subtended by any two image points on the negative plate can 
be calculated if the coordinate distances of the points can be measured and 
corrected for distortion. 

Position of Exposure Station—When the coordinate system is so assumed 
that its (X Y)-plane is tangent to the spheroid at a point of given latitude and 
longitude and the point of tangency is used as origin, the coordinates of a point 
defined by elevation, direction, and distance from the origin may be calculated 
by Eqs. 4. The straight distance between two points of known coordinates can 
be computed easily. If the space angles subtended at the exposure station by 
three points of known positions can be determined, the position of the exposure 
station which forms a tetrahedron with the three ground points is determined 
by the geometrical method or the analytical method (see subsequently in 
Example 1). 

Refraction Correction.—Since the light ray emanating from station S2 to the 
camera at station §i, Fig. 6, follows the same path as the light ray traveling in 
the opposite direction, the observed direction of point S2 from point 8; is 
deviated an amount given'by Eq. 43. To determine this error of refraction, 
the refractive constant p in Eq. 43 must be determined for the exposure station. 
As an essential requirement, at least four points, Pi, Pe, P3, and Ps—all with 
known positions—must be included in one photograph. Let X, Y, and Z be 
the coordinates of the exposure station S calculated from the space angles A12, 
Ai3, and A23 and the known positions of points P:, Pe, and P;.. From the coordi- 
nates of S, Pi, Ps, Ps, and P, and the center of curvature C of the spheroidal 
surface, the direction cosines of the lines SPi1, SP2, SP3, SPs, and SC may be 
found. The angles Ai, Az, As, and A, between the plumb line SC and the 
other four lines can be calculated and may be considered as their observed zenith 
distances, to which corrections equal to —p tan A;, —p tan As, —p tan Az, 
‘and —p tan Ax, respectively, must be added. The correction a2 for the space 
angle Ai. due to the corrections —p tan A; and —p tan Az may be found in 
terms of p by Eq. 25 on the assumption that the dihedral angle between the 

“vertical planes SCP; and SCP, is not appreciably changed by the effect of 
refraction. Similarly, the other,corrections a3, @14, G23, Ges, and sq can all be 
found in terms of p (see subsequently in Example 6). 

To account for the corrections a12, @13, and @z3, corrections 2, y, and z must 
be added to the calculated coordinates X, Y, and Z of the exposure station S 
so that Eq. 20 is satisfied. From the three equations thus obtained, the correc- 
tions x, y, and z can be found in terms of p (see subsequently in Example 7). 

Since the corrections of the direction cosines of the lines SP, SP2, and SPs 
can be found in terms of x, y, and z and therefore in terms of p by Eqs. 16, 
and ay4, de, and az4 can all be expressed in terms of p, the corrections ls, ma, 
and n, to be applied to the direction cosines of the line SP, may be found in 
terms of p by the three linear equations derived from Eq. 19d (see subsequently 


in Example 8). 


ee aie 
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The perpendicular distance D, from point Ps, which is defined by the known 
coordinates X4, Ys, and Z, to line SP, (which line, in turn, is defined by a 
point with coordinates (X + 2), (Y + y), and (2+ z) and direction cosines 
(Le +14), (Ma + m4), and (N4 + m4)) is found in terms of pby Eq.18. Because 
of the presence of other errors of measurement, it is not likely that D*, oe 0 can 
be satisfied by a real value of p; but there is always a value of p for which D*, 
is minimum. Consequently, 


may be used to determine p. Similarly if several other points Ps, Ps, and P; 
(also of known positions) appear in the same photograph, the value of p, de- 
termined from 
Abe De hee a Quinn ks i 
p 
should be more reliable. 

Equation to Determine p—Let Ls, M4, and N, be the direction cosines of 
the line which makes angles Ai, and Az with the lines SP: and SP2. When the 
coordinates of point S are changed to (X + x), (Y + y), and (Z + z) and the 
angles Ay, and Az are changed to (Ais + dis) and (Azg + ao4) to account 
for the effect of refraction, the corrections x, y, 2, l4, ms, and m4 may all be 

- found in terms of p as illustrated subsequently by Examples 7 and 8. Since 
—D4 = (X+e-—XP+ (VY ty — Yi? + Zt+e2—- 2) — [la + W) 
K(X +e — Xa) + (Ma + ms) (Y. + y — Ya) + (Na t+ 7d) (2 +2 -— Ze) PS 
is an expression involving the fourth power of p, the expression: ; 


dD*, 
Pr Oia sae cin GOA tee Ree ee (45) 
gives an equation of the form: 
atapt+azp?+as3p?=0.......... Laude oe (46), 


the solution of which gives the value of p. 
Orientation of Negative Plates—As shown in Fig. 7, point S is the rear nodal 


ait point of the camera lens; X, X’, Y, and Y’ are the fiducial marks on the nega- 


tive plate, point F is the plate intercept of a perpendicular from point S, and 
point V is the piercing point, on the negative plate, of the plumb line SV. 
Since the space angles at point S subtended by points Y, Pi, and P2 and the 
equations of the lines SP; and SP2 are determinable, the equation of line SV is’ 
computed by Eqs. 26. Similarly, the equation of SF may be found also. The 
equation of the plumb line SV is determined from the calculated coordinates of 
point S and the coordinates of the center of curvature of the spheroidal surface. — 
The angles between the three lines SY, SF, and SV are found from their _ 
equations; and the angle YFV is found from these angles by Eq. 24. Angle 
FYY’ can be found from the coordinates of points Y and F; and, as in the case 
shown, the angle between the lines FV and YY’, is equal to 


L FEY). = 7 YY! + 7 VEY, ok ene AN) 
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Since the equation of the plane of SV and SF can be determined from the 

equations of the two lines, and the equation of the plane of the meridian 

_ through point S is found from Eq. 5, the azimuth of line SF is found from the 
dihedral angle between these two planes. 

The coordinates of point F, the azimuth of line SF, and angles VSF and 
FIY’ completely define the orientation of the negative plate at the instant 
of exposure. 

As shown in Fig. 8, points S: and S» are two exposure stations and V’; and 
V’, represent the feet of the plumb lines on the spheroidal surface. The eleva- 
tions V’;S; and V’.S2 can be found from the calculated coordinates of points — 
Sand 8; by Eq. 1. The azimuth of the line SS. is found from the dihedral 
angle between the meridian plane of point §; and the vertical plane containing 
points §; and 8S». : 

Since the equations of lines S:82 and 8, V’, and the plate perpendicular §,F; 
ean be found, it is possible to compute the angle »;, Fig. 8, which may be con- 
sidered as the vertical angle of the line of sight of the camera, and the angle 
V’18i82, The horizontal angle 6: is determined from the three space angles 
at point 8; by Eq. 24. Similarly, the vertical and horizontal angles of the plate 
perpendicular S.F. at point S2 are found. These angles, together with the 
angles of orientation of the two negative plates, completely define their relative 
position in the air. 

Position of Photographed Point.—The space angles at the exposure station 
$1 made by a point P, of unknown position With known points P;, Ps, and P; 
ean be determined from their images on the negative plate by Eqs. 7 and they 
can be corrected for atmospheric refraction by Eq. 25, when the position of 
station §; and the refractive constant p are known. From these space angles 
the direction cosines of the line §,P, may be determined from Eqs. 26 or 32 
according to the number of known points available. 

Similarly, if point P, is also photographed at another exposure station S2 
and the equation of the line S2P, can be determined, the accuracy of all the 
data may be tested by computing the perpendicular distance between the two 
lines S,P, and S82P, with Eq. 21a. If this distance is reasonably small, its 
midpoint may be considered as the probable position of point P, and its coordi- 
nates may be determined by Eas. 23. 

'. Four or More Points of Known Positions in One Photograph.—In this case 
more data are available than necessary and the method of least squares is used 
to determine the coordinates of the exposure station S. 

The coordinates of point S are first calculated from the space angles of 
points P;, Pe, and P;. Then the lengths and direction cosines of lines SPi, 
‘SP., SP, and SP, are calculated; and the small corrections @12, 13, and 3 of 
the space angles Ai2, Ais, and A2s due to corrections z, y, and z which must be 
added to the calculated coordinates of point S are found from Eq. 20. 

Let the values of angles Ais, Aes, and Ass, calculated from the direction 
cosines of the four lines, be denoted by J1s, I24, and I34 and their increases in 
terms of z, y, and z, found by Eq. 20, be denoted by ‘14, t24, and ts4, respectively. 
The corrections a14, des, and a34 when added to the observed values Ais, Avs, 


ee, 2 


a 
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and Agu, respectively, must yield 


Tu + tia = Axa + RY eee PMOL MORONS te cai (48a) 

I + to4 => Aca + Bob ie Re he Size te in were eleeon Lae (48b) 
and 

I34 + 134 => Aza + gar ES ce eagle ace ela moleieaemels (48c). 
Therefore, 

ay = Ty = Au + tia SRT EA senso rele Ses ) on PRES (49a) 

aoe = Tou — Au + 124 ON Oe anes ke (496) 
and 

a34 = Is4 — Aza + 134 A Aike adh a Sty are gehen elle 1 kelfereeh ang (49c) 


which give the angular corrections d14, G24, and dss, also, as linear functions of 
x, y, and z. 

As it is possible to restrain the angular corrections if the sum Q, defined by 
the condition that— 


| Q = Wi (diz)? + Wis (a13)? + Wes (ao3)? 
+ Wis (ais)? + Woe (G24)? -+ Waa (a34)?. . (50) 


—shall be a minimum, the values of z, y, and z can be determined from the 
three normal equations derived from a =.0, a = 0, and a = 0, respectively 
(see subsequently in Example 5). In Eq. 50, W is the weight of a corresponding 
observed value. The foregoing method can be extended to other cases in- 
volving five or six ground points of known positions. 

Four or More Points of Known Positions and One or Several Points of Un- 
known Positions Common to Two Aerial Photographs.—Let the positions of four 
ground points Pi, Pe, P3, and P, be given, with that of P; undetermined, and 
_ let the available data be composed of the space angles A12, Ais, Ais, Aos, Aca, 
Asa, A1s, Ax, and Az; and A’, A’ is, A’, A’, A’, A’ 34, A's, A’ss, and A‘ 
at the exposure stations S and §’, respectively. In the solution, the coordinates 
of points S and 8’ are first computed from the given coordinates of points 
P;, Pe, and Ps by the use of the space angles at S and S’ upon these points. 
Then, corrections to these coordinates 2, y, z, x’, y’, and 2’, respectively, and 
other corrections are computed as explained in the following text. 

For the points of known positions, corrections @i2, @13, @23, @14, G24, Qs4, 
@’12, G13, Q'23, G14, Q’24, and a's, of their space angles may be expressed in terms 
of x, y, and z and 2’, y’, and 2’, respectively, by Eq. 20. For point P; (position 
unknown), the direction cosines of lines SP; and S’P; are found from their 
space angles fixed by three points of known positions such as P;, Ps, and P; by 
Eqs. 32; and their corrections I, ms, ns, l’s, m’s, and n’; are also to be determined. 
The correction ai; may be expressed in terms of 1, m1, m1, Is, ms, and ns by - 
Eq. 19; l1, mi, and n; can be obtained in terms of z, y, and z; and ns; can be | 
found in terms of J; and m; by Eq. 17. Therefore, ai; can be expressed in 
terms of z, y, z, ls, and ms. Similarly, a2, may be expressed in terms of the 
same five corrections. The correction a3; should be determined from G12, Gis, 
G23, is, and a2 by Eq. 30 and therefore may also be expressed in terms of th 

e 


_ the work of numerical calculation greatly. 


‘between photographs Gi, Ge, and G3. The method ' 
~ of successive approximation can be applied to cases 


Ha 
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foregoing five corrections. Likewise a’1s5, a’ss, and a's; may be expressed in 
terms of x’, y’, 2’, I's, and m’s. 
To make the tv lines SP; and 8’P; intersect, however, corrections 2, y, 2 


a’, y’, 2’, ls, ms, U's, and m’; must satisfy a conditional equation based on Eq. 29 


in eelich ms and n’s are replaced by their equivalents derived from Eqs. 16. 


_ By this condition of intersection, one of the corrections such as m’s can be found 


in terms of the other nine. 

The foregoing analysis shows that the basic corrections in this case are z, 
y, 2, x’, y’, 2’, ls, ms, and I’; which may be determined from a set of normal 
equations if the sum of the weighted squares of the corrections of the space 
angles is a minimum. The method is general and applicable to the case with 
(3 + qg)-points of known positions and wu points of unknown positions; but the 
number of basic corrections, being equal to 3 (2 + u) is dependent on wu only. 


- The method can also be modified to include such data as A4; and A’4;. If the 
refractive constant p is known, the space angles should first be corrected for 
_ refraction by Eq. 25, and then they should be used as the data to start the 


aforementioned method of analysis. 

If (because of some special arrangement to be explained subsequently) the 
position of the exposure station S has been determined previously, the correc- 
tions x, y, and 2 may be considered equal to zero. The number of the basic 
corrections is then decreased by 3, thus decreasing 


Analysis of Partly Overlapping Aerial Survey.— 
Fig. 9 represents an aerial survey with a 60% overlap 


of this kind. It is assumed at the beginning that: Fra. 9 
(1) At least three points of known positions, such as 
Pi, Ps, and Py, are available for the determination of the exposure station 8; at 


- which photograph G; is taken; (2) in photograph Gz only two points of known 


positions such as P; and P, can be identified; and (3) each of the succeeding 


- photographs contains at least two well-defined points in common with two 
_ previous photographs. It is also assumed that: (4) The angles of tilt are small 
and (5) the altimeter records are reasonably reliable. Under these assumed 
conditions, the exposure station may be represented approximately, in plan, 


by the principal point, and in elevation by the altimeter reading. If the prin- 
cipal point of a vertical photograph can be identified in a previous photograph 


_ by its surrounding topographic details, its position in plan is estimated or de- 


termined by the radial-line method of plotting. 

From the known coordinates of points P3 and P., and from the estimated 
coordinates of the exposure station S2 at which photograph Gz is taken, the 
direction cosines of the lines S2P3 and S2P, and their including angle J’’3, may 
be calculated. The change 7’’34 of I’’34 due to corrections 22, yz, and 22 to be 
added to the estimated coordinates X2, Ys, and Zz of point S2 may be found by 
Eq. 20. Since the sum (I's, + 7’’s4) must be equal to the observed value A’’ss 


of the space angle at point So, 


18 AERIAL SURVEY 


which is a linear equation in 22, yz, and 22. If it is believed that z2 and y2 are 
relatively small because X2 and Ye are more reliably determined than Zz, 
Eq. 51 can be used to find z2 by omitting the terms involving z2 and y2, In 
that case, the calculation should be begun again with X2, Y2, and (Z2 + 22) 
as the initial values. It is essential that good initial values be obtained by a 
simple method so that the entire process of calculation can be shortened. 

The direction cosines of line S2P; may be calculated from the computed 
values of the direction cosines of lines S2P3 and S.P, and the space angles A’’ss 
and A's; by Eqs. 26. The changes in the direction cosines of lines S2P3 and 
S.P, due to the corrections x2, y2, and 22, and the corrections I’’s, m’’s, and n’’s 
to be applied, to the calculated values of the direction cosines of line S2P; must 
not alter the values of the space angles A’’3; and A’’4;. Therefore, 


L's ’’, a Loe U's ita Mir m’’s ae M"; ms +N"; n''s as Ns n''s a0. (52a). 
; L* i’, af Li’, ee aa M"; ms a M", ms = N"; n''s a N"; ns aside .(52b) 


and 
Ls U's + M". ms + N"s ns Sa ies che teh ee ate (52c) 


in which all the capital letters denote numerical constants determined from 
known or assumed values. Corrections l’’3, m’’3, 1'’3, 1’’4, m’’s, and n’’,4 can 
be found in terms of 22, y2, and z2 by Eqs. 16; therefore corrections l1’’s, m’’s, 
and n”’; can also be obtained in terms of 22, y2, and 2. by Eqs. 16. The line 
S.P; must intersect line 8:P;, however; and the equation of the latter has been 
determined completely in the analysis of the previous survey. Hence, the 
equation of intersection based on Eq. 22 involves only 22, yo, 22, l’’s, m’’s, and 
n'’; as the unknowns; and, as the latter three corrections may be found in terms 


of the corrections of coordinates, the final result is again a linear equation in — 
“2, Y2, and z:. Similarly, the intersection of the lines S,:P. and S2P¢ gives an- ~ 
other equation in 22, y2, and 22. These three equations are used to determine » 


the values of the three corrections. 


As the errors of the initial values Xe, Ye, and Z. may be very large, the . 
corrected values (X2 + 2), (Y2 + yo), and (Z_ + 22) are usually only better — 


approximations. With these revised values of the coordinates of point Se as 
the new base for starting, the foregoing procedure is repeated (see subsequently 
in Example 9). 

Since there may be other well-defined points such as P; common to photo- 
graphs G, and Gz, and since from each common point a linear equation in 22, 


y2, and 2: can be derived, more equations may be obtained than are required — 


for the determination of the three unknowns. In such a‘case, in the second 
course of approximation, the method of least squares is used. 

After the coordinates of point S: are thus calculated, the equations of the 
lines from point S, to all points in photograph Gz may be computed, and the 


positions of the points common to photographs G, and G» may be determined 
by solving Eqs. 23. 


The same method of calculation may be used for obtaining the positions of — 


the succeeding exposure stations and points, until the survey ends or closes on 
another point of known position. 
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Aerial Triangulation —Present practice in aerial surveying is mainly based 
on the use of partly overlapping (generally 60%) exposures as shown in Fig. 9. 
The method of computing such a survey is not straightforward although results _. 
may be obtained which are as reliable as the angular data. 


DXDT 
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A fully duplicated, and also partly overlapping, aerial survey may be made 
with two cameras so arranged that their fields of view have a longitudinal 
overlap of 20% as shown in Fig. 10(a) or of 60% as shown in Fig. 10(6). When 
cameras arranged in such a combination are operated simultaneously by electric 
control so that the exposures of each camera have a 60% overlap of their own, 
each photograph of one camera will be nearly duplicated by some photograph 
‘of another. 

Let a series of triangles Ti, To, ---, T, on a strip of ground surface (see 
Fig. 10(c)) be so spaced that;adjacent pairs will appear together in the same 


TABLE 14.—ScHEDULE oF SE- TABLE 1B.—ScuHEDULE oF SE- 
QUENCE IN THE Use oF TRI- QUENCE IN THE Use oF TRI- 
ANGLES IN AERIAL TRIANGU- ANGLES IN AERIAL TRIANGU- 

LATION COMPUTATIONS; LATION COMPUTATIONS; 
25% OVERLAP 60% OvERLAP 
Exposures Backsights ; Foresights Exposures Backsights Foresights 
Ei Ti Te Fi Ti Te 
EH’: Ti Te E's Ti Te 

Ts Wy Ee Te Ts 
Be Ts T; HE’: Tz Ts 
E oT: That E: Ta Tas 
B’ny2 | T. | Bott Bout Ta Th41 


photograph; and let the simultaneous exposures of the two cameras at the jth 
exposure station be denoted by E; and H’;. If triangle T, is composed of points 
of known positions and may be used to determine the position of the exposure 
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station, the two sets of exposures can be compiled as in Tables 1A and 1B, , 
showing the sequence in which calculations are made. As in the case of a, 
ground survey, the sights used to determine the positions of the instrument are | 
called the ‘‘backsights,” and those used to determine the points sighted are : 
called the ‘‘foresights.” In either case the backsights of triangle T, are used . 
to find the positions of the exposure stations 8; and 8’; or 8; and S’2 and the : 
foresights from two stations of known positions are used to find the vertexes | 
of triangle T2. This process of calculation is continued until the series of ' 
triangles ends on another point of known position. 

As the relative position of the exposure stations S, and S’, is identical to | 
that of the two cameras in the airplane, when their positions are determined 
independently, additional checks on the accuracy of the data and the numerical 
calculations are provided by the spacing of the cameras. 

Error of Closure of Aerial Triangulation—When an aerial survey begins 
from a ground triangle composed of points of known positions and ends on 
another known ground point, the coordinates Xa, Ya, and Z, of the terminal 
point determined by the aerial survey should agree with the more reliable 
values X,, Y,, and Z, determined by the ground survey. The errors in three 
directions are: er = Xa — Xq, €y = Ya — Yq, and e. = Za — Zz, respectively, 
and the error of closure is equal to e, = Vez + e, +. Let the distance 
between one of the starting points and the terminal point be D. The ratio 


5 , reduced to the form 5 , will serve as a measure of the accuracy of the survey. 


If it is assumed that the errors e., ey, and e, are accumulated equally from 
all parts of the survey, the corrections to be added to the calculated coordinates 


of an intermediate point should be: meee : ae and — pe , Tespec- 


tively, in which D; is the dis- 
tance of the point from the 
starting point. 


TABLE 2.—CoorpI1naTEes oF GROUND 
Points, IN METERS 
(1 Meter = 3.28083 Ft) 


ILLUSTRATIVE EXAMPLES 


Point x Y Z 
Boo 2578 re 3.876 88 ese satin 2° ae in Table : 
Tae 769: 643.21 135.27 ar 
Pats cee 6,043.75 5,872.15 146.73 heey cir te ube inal gris) 
Pe: cadd. 2/385.00 5,934.00 76.30 making sample calculations that 


——— —————————— _ will illustrate the general fea- 

tures of the methods and offer 
some indication as to the accuracy of the results to be expected. The positions 
of two exposure stations are also assumed so that the field angles to the various 
ground points are, for station Si, 


Angle Value 
BiB Pas cee eh a cocea rh ea 67° 31’ 
Piss Roi eee hh On ee ee 69° 47’ 
far ot fauea ee © AGN oe MR ie Ay > 75° 38’ 
BAP 4.4 vat See bis, es pon ce eee 29° 08.011’ 
Pasi Page od ke scien oe eka Dee 89° 53.239’ 
Ras Pa cusi Uvi hit OG Onan eee 55° 47.648’ 


ar 
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and, for station S.2, 


Angle Value 
denool rete RAC er Sidr eas Sst: 44° 04.508’ 
i age Vath et So oa ORS ey Ae ea a soe ae 62° 05.500’ 
1 ESSN SAN pec 0 gle A a a 29° 24.033’ 
Mookie enya MR Taian tos he lcs con he oe Ss 49° 10.383 
AERIS, TEA BSS oa nab ce ene a ee ee a 68° 52.333’ 


Example1. To Find the Coordinates X1, Y1, and Z; of Exposure Station 8;.— 
The lengths D1, Do, and D; of lines 8;P;, 8:P.2, and 8,P; must first be determined 
from the equations: 


D*, — 2(0.8824147) Di De + DY, — 15,173,542 = 0....... (5a) 

D*%, — 2(0.3455712) D: Ds + D%; — 15,997,524 =0....... (536) 
and 

D*, — 2(0.2481263) D2 Ds + D%, — 17,959,339 = 0....... (580) 


The results are summarized in Table 3. The final values are accurate to the 


- TABLE 3—Computation oF Lenetus D sy Eas. 53 (IN MzTERs) 


No. Assumed | Corrected | Computed | Corrected | Computed | Corrected | Computed |Computed 
(1) (2) (3) (4) (5) (6) (7) (8) 
ae Me +3,544.60 | +3,504.87 | +3,465.27 | +3,487.09 | +3,447.09 | +3,424.18 | +39.6 


545 +3,543.957| +-3,463.818] +3,464.918] +3,445.691| +3,446.791) +3,466.005} — 1.043 
“| $3" ,543.957| +3,543.963| +3,464.916] +3,464.910| +3,446.796| +3,446.791| +3,464.904| + 0.0057 


a 


third decimal place. The coordinates of point 8; (in meters) are then found 
to be: Xi = 4,953.549, Yi = 3,827.388, and Z; = 2,698.353. 

Example 2. To Find the Elevation of Station S:.—Let the radius of curvature 
of the spheroidal surface be R = 6,363,200 meters (1 meter = 3.28083 ft). 
It is obvious that the third term 


in Eq. 1 can be neglected in this TABLE 4.—CompuTaTION OF 
case and that the elevation P;Si, DirEcTION COSINES 
; 19,593,000 
therefore, is 2,698.353 + <3 5 
4 fe + 6 365,900 __Line L M N 
= 2701 SiP;....| —0.0007781 | —0.0006012 | —0.9999995 
2101.43) mi. 4 E 5 SiPi Sy. —0.6701815 | +0.0138241 | —0.7420682 
Example 3. To Find the Direction - SiP2..... +0.2354551 | —0.6303707 | —0.7397255 
: : SiPiccn. 40.3162945 | +0.5932365 | —0.7402895 
Cosines of Lines SiP;, SiP1, SiP2,  SiPals::: —0.6069250 | +0.4977735 | —0.6195673 


3:P3, S:iP4, and the Included Angles.— 
The direction cosines of the plumb 
ine S,P; are determined from the coordinates of station S:, and the center of 
survature (X = 0, Y =0, and Z = — 6,363,200). Similarly, the direction 
osines of the other four lines may be found from the known coordinates (see 
Table 4). From the direction cosines the following space angles may be com- 
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puted by Eq. 19a: 


Angle Value 
PAP ace Se a Pee ee es eee 42° 02.893’ 
PSPs aa ae aap eens 42° 16.518’ 
PSPs eRe awn ae ates Seles Py 42° 17.720’ 
BSjP eo nies ee ek ees oe ee eee ph 42.172’ 

and, as measured to Pu, 

Angle Value 
RR PAU yi uh ee elem hae ee 29° 08.681’ 
RE Pip ee. alyssa Nh ioe le te 6 ch aden okay De ot 89° 54.417’ 
Be oy APPR Se RM Re RR Pag a Shiboy Pager roe? 55° 48.389’ 


Example 4. To Find the Direction Cosines of Line SiP4 from the Direction 
Cosines of Lines 8:P1, S:iP2, and SiP3—The direction cosines of line $:P, found 
from the calculated values of angles 

TABLE 5.—D1rrerENces BETWEEN’ P,S,P,, P2SiPs, and P3:S:P, by Eqs. 
THE ASSUMED AND COMPUTED 32 are: Ly = — 0.6069253, M, = + 
VALUES OF ANGLES AT 0.4977740, and Ns = — 0.6195678. 
Sration §; ' Example 5. To Find the Most' 


————— | Probable Position of S:.—-Tableme 

Angle Assumed | Computed Error demonstrates the differences between: 
PiSiPs....| 29°08.011’ | 29° 08.681’ | 0°00.670' the assumed and the calculated values 
PSiPi....| 89° 53.239’ | 89°54.417’ | 0° 01.178" 


PiSiPs.. | .| 55° 47.648’ | 55°48.389’ | 0°00.741' of the angles PiSiPs, P2SiP.s, and 

———_ PS, Py. If these errors are to be! 

distributed to the six space angles at: 

station §, by introducing corrections 21, y:, and z; to the calculated coordinates 
of station S,, it may be found from Eq. 20 that 

a2 = — 0.30007 x; — 0.89364 y: — 0.97174 2.......... (54a)' 


a3 = — 0.26216 x, + 0.40511 y: — 1.0170 2;........... (546) 
Qo3 = 0.42573 i 0.03279 a 1.1427 3s ee a a OE Re (54c) 
and Egs. 20 and 49 will yield 


and 


di4 = 0.670 — 0.27669 x1 + 0.26535 yi — 0.27829 2, 
de, = 1.178 — 0.41046 x, — 0.01784 y; — 1.2186 z,........ (55d). 


0.741 — 0.30059 2, + 0.50608 y: — 0.63048 z....... (55¢) 


and 


34 


The angular corrections are all reduced to minutes of are by multiplying the: 
results by 3,437.7, which is the number of minutes in one radian. The normal 
_ equations derived from these corrections, with equal weights, are 


0.6786 x, — 0.2203 y; + 0.8364 z; = 
— 0.2203 2; + 0.6470 y, — 0.3633 2, 


0.8364 21 — 0.3633 y: + 5.2323 z, 


| eae | 
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from which z; = 0.9105, y: = — 0.3851, and z, = 0.2258. These are the cor- 
rections to be added algebraically to the coordinates of Sj, respectively. 

Example 6.’ To Find the Corrections of the Space Angles for Atmospheric 
Refraction.—By using the calculated values of the space angles between the lines 
SiP;, SiP,, SiP., S:iP3, and SiP,, Eq. 25 yields: ai2 =— 1.3368 P, &3 = — 1.3948 P, 
G23 = — 1.5523 p, ais = — 0.5873 p, do4 = ~ 2.0283 p, and a3, = — 1.0970 p— 
in which p is in minutes. 

Ezxample?. To Find the Corrections of the Coordinates of 8: in Terms of p.—By 
equating the errors of refraction in angles P,;SP2, P:SiP3, and P.S,P3 to those 
due to the corrections 21, y1, and 2; to be added to the calculated coordinates of 
station S,, three equations (in meters)— 


2 0.30007 x: + 0.39364 y: + 0.97174 2: = 1.8368 p....... (57a) 
: 4 0.26216 2, — 0.40511 y: + 1.0170 z, = 1.3945 p........(57b) 
an 

— 0.42573 21 + 0.03279 yr + 1.1427 2 = 1.5523 p........(57c) 


—may be obtained which give: z; = 0.0231 p, y1 = 0.0042 p, and z; = 1.3669 p. 

Example 8. To Find the Corrections of the Direction Cosines in Terms of p— 
From the values of 21, yi, and 41 
found in Example 7 and from 
other known values, by Eqs. 16, 
the corrections of the direction 


TABLE 6.—CorRREcTIONS FoR DiREC- 
TION CosINES, EXAMPLE, 8 


: Li l 
Bes Ol CaP yr his sat Or § yt ee 
may be determined as shown in Biber ae po007e2 > $o.0001822 p =3.00 0001703 p 
T : : ak: —0.0000985 p | —0.0001737 p 13 p 
Table 6 in which U4, ms, and ma gipirst:: +40.0001706 p | —0.0001457 p | —0.0002847 p 


are found from the three linear 
equations derived from Eq. 19b. 

_ Example 9. To Find the Coordinates of Station S2 from Space’ Angles, and 
the Known Positions of Points P; and P:—The initial values assumed for the 
coordinates of station S2 are: X2 = 4,400, Y2 = 6,210, and Z. = 3,160. With 


- TABLE 7.—Coorpinates oF Station S2 COMPUTED FROM SPACE ANGLES 


Xe Ye Z2 
ae As- | Correc-| G tea | AS, | Correc-] o ted | _As-, | Correc-| Corrected 
% sumed | tions orrected | sumed | tions orrected | sumed | tions 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 
1 Isao +4,400 | —58.5 | +4,341.5 |+6,210 | —69.7 | +6,140.3 |+3,160 | —19.8:| +3,140.2 
J 2.....)+4,341 | +17.6 | +4,358.6 |+6,140 | +10.4 | +6,150.4 |+3,140 | — 8.2 | +3,131.8 
True value..| .... «se. | +4,358.00) .... -.+. | +6,150.00| .... ...+ | +8,182.00 


these and other known values the solution of Eq. 51 yields: 536.9 x2 — 2,160.7 ye 
— 1,499.2 2. = 150,214, as a correction for angle P,S2P2. Similarly, Eq. 22 
yields 
0.1691 x2 + 0.1015 y2 + 0.5552 22 = — 28.02 
znd 

0.2796 x2 — 0.5559 y2 — 0.4415 22 = 31.12. 


ae 
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for making lines S2P3 and S2P, intersect the lines $:P3 and S,P4, respectively. . 
The solution of these three equations gives: 42 = — 58.5, y2 = — 69.7, and 
z. = —19.8. The foregoing results and those of a further process of approxi- 
mation are shown in Table 7 which indicates a steady convergence to the 
true values. 
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